Comment on fractality of quantum mechanical energy spectra 
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The fractal properties of the energy spectra of quantum systems are discussed in connection with 
the paper by Saiz and Martinez [Phys. Rev. E 54, 2431 (1996)]. It is shown that for discrete energy 
levels the Hausdorff-Besicovitch dimension is zero and differs from the Renyi scaling exponents 
computed by the standard box counting algorithm. The Renyi exponents for the inverse power 
series data sets (x n = n = 1,2,...) are computed analytically and they are shown to be 
— and, as a consequence, do = 1/3 for the Balmer formula. 
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In last years there is an increasing interest in searching 
for fractal structures in physics and in looking for frac- 
tal signatures of chaos at the quantum mechanical (QM) 
level [HHI), as well. In fact, a simple QM system with 
fractal energy spectrum has been found long time ago 
H . Similar structures have been found in recent years in 
other systems that are of great practical importance: the 
quasiperiodic semiconductor microstructures, the quan- 
tum Hall effect and the Anderson localization In 
this paper I would like to comment on possible fractality 
of QM spectra and I will discuss Renyi exponents [|]|]] 
for some special types of point spectra generated by sim- 
ple probability distributions, as well as for Balmer like 
energy levels computed in M 



n > m 



1,2,3,4,5. (1) 



Typical QM systems have either discrete point spectra 
(localized states, like for the harmonic oscillator), contin- 
uous spectra (extended states, like for the free particle) or 
both (localized and extended states above some thresh- 
old energy, like for the hydrogen atom). In addition, the 
models have been found that have neither extended nor 
localized states and their energy spectra have been shown 
to be fractal These are models that describe infi- 

nite crystals (periodic or almost periodic). In this paper 
I will limit discussion to the case of discrete spectra. 

For the discrete spectra it has been shown that their 
nearest neighbor spacing (NNS) of energy levels has Pois- 
son or Wigner probability distribution (for correspond- 
ing chaotic Hamiltonians) jl0],[ll]]. Such spectra cannot 
lead to any fractal structure as any reasonable probabil- 
ity distribution cannot give the Cantor like structures (in 
particular, the Poisson, Wigner or Thomas-Porter dis- 
tributions are regular enough). On the other hand, the 
energy spectra with finite number of accumulation points 



cannot be generated by a reasonable (regular) probabil- 
ity distribution. However, using the correct definition of 
the fractal (Hausdorff-Besicovitch) dimension |l^,|l3| one 
can easily show that for any set with finite number of ac- 
cumulation points the Hausdorff-Besicovitch dimension 
is zero, dn = 0. Hence, for any nuclear or molecular dis- 
crete energy spectra we should get their fractal dimension 
equal to zero. 

In practical (numerical) calculations of the fractal di- 
mension the box counting algorithm and the Renyi ex- 
ponents (improperly called "dimensions") are used ^4 
The whole interval is divided into N equal subinterva 
("boxes") with rii(N) (i — 1,2, ... ,N) data points in 
each box. Defining the measure Pi(N) = rii(N)/ntot, 
where n to t = Si n «(^) ^ s the total number of data 
points, we have the following definition of the Renyi ex- 
ponents 



1 Hm InEiPfW 



q—lN- 



hxN 



(2) 



Here, for q — (the capacity "dimension") we assume 
summation over non -empty boxes only, and (^) becomes 
equivalent to 



In M (A) 
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(3) 



where M(N) is the total number of non-empty boxes. 
For q = 1 (the information "dimension") the de l'Hospital 
rule is applied to the formula (0) and one gets 



d\ = — lim 

AT->-o« 
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(4) 
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Once again, it is important to remember that the Renyi 
exponents are not the (fractal) dimensions and for some 
"pathological" sets they may differ from the Hausdorff- 
Besicovitch dimension [ |12|Jl^ 1. In particular, this hap- 
pens for the inverse power like discrete sets, as will be 
shown below. However, even for such pathological sets 
one can use the Renyi exponents just to describe the scal- 
ing properties of the data. With this remark in mind we 



1 



prove the following 

Theorem: For the inverse power series x n = l/n a , where 
n = 1,2,3,... and a > the Renyi exponent do is 



do 



1 



1 



(5) 



Proof: Let us denote by n sng i the number of points 
that are single in their boxes. We have M(N) > n sng i 
and, from the condition x n — x n+ \ > 1/N (distance be- 
tween such points is greater than the box size = 1/-/V), 
we get n sng i — a T +° : N T +^. In fact, in the last equa- 
tion there should be taken integer part of the right hand 
side but this is unimportant as we are interested in the 
limit N — > oo. Using the definition ||with q = one gets 
asymptotically for N — > oo the lower limit do > (a+l) -1 . 
To get the lower limit we took into account points in the 



UK 



sngl 



To obtain an 



interval [x sn gi, 1], x sng i — %n sn 

upper limit we have to take into account the remaining 
boxes as well. First, let us notice that the first box, i.e. 
the interval [0, 1/N] contains infinite number of points 
for any N, as zero is the accumulation point, and the in- 
terval [1/N,x sng i] contains at most (x sng i — 1/N)/(1/N) 
non-empty boxes. This gives the following upper limit 
for M(N) 



M(N) < a~ 
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and, as a consequence, we have 



In 
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a — (i + I) n— 
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Hence, we have proven the eq. (g). 

This result gives us the scaling exponent do = 1/2 for 
the harmonic series (as was stated in while for the 
Balmer like series (a = 2) we have do = 1/3 ~ 0.33, in 
contrast to the result presented in ||, where the value 
0.61 was obtained numerically. Of course, neither the 
Hausdorff-Besicovitch dimension nor the Reny exponent 
be can changed by the overall shift of the whole set (due 
to the 1/m 2 term in (Q)), or by rejection of any finite 
number of the data points (n > m), or even by super- 
position of any finite number of such sets (as we have 
in (g) logarithm in the numerator). In fact, superposi- 
tion of two fractal sets has the fractal dimension equal to 
the maximum fractal dimension of the constituent sets 
|l3f . Hence, our result cannot be changed by any simple 
composition of Balmer like sets. 

I have done some numerical checks of the analytical re- 
sult ([5]) with the box counting algorithm and both meth- 
ods are consistent. In particular, for a — 2 and 1000 
data points the box counting method gives: do — 0.321, 
where 6 points fit to the straight line with x 2 ~ 0.0072 
and the correlation parameter r ~ 0.99988 (the finest 
division was N — 2 25 ). 

Finally, I would like to mention that the modified scal- 
ing exponents were suggested to investigate QM systems 



with discrete spectra [^Jj]. Namely, the probabilities 
Pi{N) obtained by the simple level-counting in a box 
have been changed by summation of the de-excitation 
probabilities. In this case, not only the information con- 
tained in the energy spectrum but also in eigcnfunctions 
is employed. However, even in this case, up to now, we 
are unable to find any interesting physical interpretation 
of the scaling exponents. 
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